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ABSTRACT 

We use images of high spatial and temporal resolution, obtained with the Rapid Oscillations in the 
Solar Atmosphere instrument at the Dunn Solar Telescope, to reveal how the generation of transverse 
waves in Type I spicules is a direct result of longitudinal oscillations occurring in the photosphere. 
Here we show how pressure oscillations, with periodicities in the range 130 - 440 s, manifest in small- 
scale photospheric magnetic bright points, and generate kink waves in the Sun's outer atmosphere with 
transverse velocities approaching the local sound speed. Through comparison of our observations with 
advanced two-dimensional magneto-hydrodynamic simulations, we provide evidence for how magneto- 
acoustic oscillations, generated at the solar surface, funnel upwards along Type I spicule structures, 
before undergoing longitudinal-to-transverse mode conversion into waves at twice the initial driving 
frequency. The resulting kink modes are visible in chromospheric plasma, with periodicities of 65 - 
220 s, and amplitudes often exceeding 400 km. A sausage mode oscillation also arises as a consequence 
of the photospheric driver, which is visible in both simulated and observational time series. We 
conclude that the mode conversion and period modification is a direct consequence of the 90 degree 
phase shift encompassing opposite sides of the photospheric driver. The chromospheric energy flux of 
these waves are estimated to be «3xl0^ Wm^^, which indicates that they are sufficiently energetic 
to accelerate the solar wind and heat the localized corona to its multi-million degree temperatures. 

Subject headings: methods: numerical — magnetohydrodynamics (MHD) — Sun: atmosphere — Sun: 
chromosphere — Sun: oscillations — Sun: photosphere 



1. INTRODUCTION 

The origin of the energy required to power the solar 
wind, and heat the Sun's corona to its multi-million 
degree temperatures, remains an elusive mystery in 
physics. One promising mechanism involves the cre- 
ation of Alfvenic waves near the solar surface, which 
can penetrate upwards through the Sun's atmosphere 
with minimal reflection or energy loss ()Erdelvi fc Feduiil 
l2007f) . To date, there has been great controversy sur- 
rounding the interpretation of periodic transverse mo- 
tions as Alfven waves. Magneto-hydrodynamic (MHD) 
waves with a wave number, m — 1, can be consid- 
ered "kink" oscillations when structured by a cylindri- 
cal waveguide. The associated periodic transverse mo- 
tions imply a degree of plasma non-uniformity, result- 
ing in intrinsic MHD waves exhibiting mixed (i.e. slow, 
fast, Alfven) properties depending on the local plasma 
parameters. As a result, the term "Alfvenic" was intro- 
duced to describe waves which are predominantly charac- 
terised by the signatures displayed by pure Alfven waves. 

|d.jess@qub.ac.uk| 



Accordingly, transverse kink oscillations may often be 
considered Alfvenic in nature providing they are in- 
compressible, exhibit no intensity fluctuations along the 
structure, and display periodic displace ments with the 
magnetic tension as the restoring force (|Goossens et al.l 
[2009). Thus, by deflnition, Alfvenic waves require strong 
magnetic field concentrations, and/or a steep density gra- 
dient with the external plas ma, to act as waveguides 
(jVan Doorsselaere et al.l[2008l ). allowing energy to be di- 
rectly channelled through the solar atmosphere. In the 
Sun's atmosphere, magnetic field lines clump together 
into tight bundles, forming flux tubes. The ubiquitous 
nature of Alfvenic waves in magnetic flux tubes has been 
demonstr ated in a range of chromospheric and coronal 
plas mas (lOe Pontieu et al I l2007bl: iTomczvk et aH 120071 : 
Baner iee et al.ll2009t iJess et al.ll2009[ r 

Spicules are dynamic, straw-like magnetic structures 
found in the solar chro mosphere, and can be d ivided 
into two distinct classes ([Zaaarashvili fc Erdelvil [2009'). 
We will focus on Type I spicules, which are ubiq- 
uitous throughout the solar atmosphere, and longer 
lived (ftilO min) than their je t-like Type II counter- 
parts (|De Pontieu et al.l l2007a|) . The narrow shape of 
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Fig. 1. — Simultaneous images of the Ha core (chromosphere; upper left) and G-band (photosphere; upper right), acquired at 14:06:21 UT 
on 28 May 2009. Dashed white lines in the Ha core image highlight the spatial positions where time-distance cuts are made. The blue 
trace in the G-band image denotes the position of a chromospheric spicule (visible in the left panel as a dark, straw-like structure), where 
one end is anchored into the photosphere above a MBP. A red contour indicates the location where a high concentration of longitudinal 
oscillatory power is present. The lower panel is a sample Ha time-distance cut, obtained 4000 km (5.5 arcsec) from the underlying MBP, 
revealing an abundance of periodic transverse motions in the solar chromosphere (see also Supplementary Movies 1 & 2). White arrows 
highlight a trough and a peak of a typical transverse oscillation. The scale is in heliocentric coordinates, where 1 arcsec approximately 
equals 725 km. 



all spicule structures (<600 km), coupled with their 
rapidly evolving characteristics, places an utmost re- 
quirement for both high-spati al and -temporal r es olution 
inst rumentation. Recently, IDe Pontieu et al.l ()2007bO 
and IHe et al.l (|2009) observed Type I spicules, which dis- 
played periodic transverse motions with periods ranging 
from 50 - 500 s, coupled with velocity amplitudes ex- 
ceeding 25 kms""'^. Unfortunately, such image sequences 
of spicules have been tied to limb observations. While 
observations at the solar limb allow spicule height in- 
formation to be easily retrieved, they do not allow si- 
multaneous measurements of the Sun's photosphere to 
be studied due to inhere nt li ne-of-sight effects. How ever, 
IDe Pontieu et al.l (|2011f ) and lMcIntosh et al.l (|20lH) have 
recently shown how spicule dynamics make them ideal 
candidates in the quest for discovering efficient energy 
propagation into higher layers of the Sun's atmosphere. 

While transverse kink motion has been shown to 
be abundant in spicule observations, the underlying 
cause of the periodic motions has remained specula- 
tive. Overshooting of convective motions in the photo- 



sphere, granular buffeting, rebound shocks, and global 
p-mode oscillations have all been suggested a s candi- 
dates for the creation o f spicul e oscillations (IRobertsI 
[1971 [Sterling fc HoU^g 1988; Vranies et all 120080 . It 
is known that oscillations in the lower solar atmosphere 
are prone to aspects of mode coupling, whereby one 
particular magneto- acou stic wave mode (normally clas- 
sified as 'fast' or 'slow'; iNakariakov fc Verwichtel I2005D 
may couple, and hence transfer wave energy, into an- 
other mod e when a tmospheric conditions are suitable 
(UlnischncidereLal. 1991; Kalkofcn 1997). The rapid 
drop in density above the solar surface, coupled with the 
close proximity to oscillatory drivers, makes the photo- 
sphere and chromosphere an ideal locatio n for the occur- 
rence of efficient mode coupling. Indeed, iMcAteer et alj 
■: ■_*()();'); uncovered upward-propagating transverse oscil- 
lations, generated in large (>2000 km diameter) pho- 
tospheric bright points, which coupled into longitudi- 
nal waves at chromospheric heights. Often, the cou- 
pled oscillations demonstrate power at frequencies equal 
to twice the original, as a result of non-linear interac- 
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TABLE 1 

Spicules observed in both Ca ii K and Ho datasets. 
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Velocity 


MBP Anchor 


Number 


(km) 


(s) 


(km) 


(kms-i) 
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670 
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Spicule used for display purposes in Figures [TlfSl 
Spicule used for display purposes in Figure l4l 



tions (jUlmsclineider et al.lll99li ). Contrarily, it has been 
demonstrated, both analytically and theoretically, that 
when the magnetic pressure is approximately equal to 
the gas pressure (plasma /3 = 1) in a solar flux tube, 
longitudinal to transver se mode coupling may also occur 
pe Moortel et al.|[200l ). To date, however, this form of 
"reversed mode coupling" has not been verified observa- 
tionally. In this Letter, we undertake a multi-wavelength 
approach, allowing Type I spicules to be traced back to 
their anchor points on the solar surface. We provide 
conclusive evidence on how these ubiquitous oscillations 
manifest, and suggest what role they may take in the 
localized heating of the solar corona. 

2. OBSERVATIONS 

The data presented here are part of an observing 
sequence obtained during 13:46 - 14:47 UT on 2009 
May 28, with the Richard B. Dunn Solar Telescope 
(DST) at Sacramento Peak, New Mexico. The newly- 
commissioned Rapid Oscillation s in the Solar Atmo- 
sphere (ROSA; IJess et aI]|2010bD six-camera system was 
employed to image a 69.3" x 69.1" region positioned at 
solar disk center. A spatial sampling of 0.069" per pixel 
was used for the ROSA cameras, to match the telescope's 
diffraction-limited resolution in the blue continuum to 
that of the CCD. This results in images obtained at 
longer wavelengths being slightly oversampled. However, 
this was deemed desirable to keep the dimensions of the 
field-of-view the same for all ROSA cameras. 

During the ob servations presented here, high-order 
adaptive optics (|Rimmelel |200j) were used to correct 
wavefront deformations in real-time. The acquired im- 
ages were further improved through speckle reconstruc- 
tion algorithms (jWoger et al.l 120081 ) . utilizing 16 1 
restorations, providing reconstructed cadences of 0.528 s 
(G-band) and 4.224 s (Ca ii K and Ha). M ore details 
of our instrumentation setup can be found in iJess et all 
(|2010aD . Atmospheric seeing conditions remained excel- 
lent throughout the time series. However, to insure ac- 
curate coalignment in all bandpasses, broadband time 
series were Fourier co-registered and de-stretched using 
a 40 X 40 grid, equating to a ^1-7" s eparation between 
spatial samples ()Jess et al.ll2007l I2008D . Narrowband im- 
ages were corrected by applying destretching vectors es- 
tablished from simult aneous broadband reference images 
(|Reardon et al.ll2008[) . 



3. ANALYSIS AND DISCUSSION 

The ROSA field-of-view shows a range of features, in- 
cluding exploding granules and a multitude of magnetic 
bright points (MBPs; Fig.[T]). A large conglomeration of 
MBPs were located at heliocentric coordinates (4 arcsec, 
—3 arcsec), or SOI WOO in the solar north-south-east-west 
coordinate system, providing an opportunity to exam- 
ine these kiloGauss structures without the line-of-sight 
spatial offsets normally associated with near-limb obser- 
vations. Examination of Ha images reveals a wealth of 
chromospheric jets and spicules, anchored into the solar 
surface at the locations of MBPs. These data reveal how 
un-anchored ends of magnetic flux tubes are observed 
to sway periodically over many complete cycles (Supple- 
mentary Movies 1 & 2). To quantify the oscillatory pa- 
rameters, a series of one-dimensional slits, each separated 
by 1 arcsec (725 km), were placed perpendicular to the 
Type I spicules along their entire length (upper panels 
of Fig. [T]). The resulting time-distance cuts reveal nu- 
merous periodic transversal displacements of these chro- 
mospheric structure s as a function of ti me. The sinusoid 
fitting algorithms of lBallai et al.l (|2011D were used to ex- 
tract periodic signatures from the transversal displace- 
ments present in the time-distance plots. We concentrate 
on the seven most prominent Type I spicules, observed 
in both Ca ii K and Ha time-series, which can be ac- 
curately traced back to their photospheric counterparts. 
The lower panel of Figure [T] displays a sample time- 
distance cut, revealing how these spicules demonstrate 
transverse kink oscillations above their corresponding 
MBP in the period range 65 - 220 s, with displace- 
ment amplitudes reaching as high as 670 km. Veloc- 
ity amplitudes often exceed 15 kms~^, indicating these 
dynam ic motions are close to th e chromospheric sound 
speed ()Fossum fc CarlssonI I2005D . Table [1] displays key 
characteristics of the spicules under investigation. 

No periodic transverse motions in either the MBP 
edge, or bright central peak, were detected in simul- 
taneous photospheric G-band images, implying that ei- 
ther spicule motion is generated independent of processes 
manifesting on the solar surface, or that a form of effi- 
cient mode conversion between the two layers (photo- 
sphere and chromosphere) is occurring. To investigate 
the latter, G-band intensities were averaged over the en- 
tire MBP structure under investigation, with fast Fourier 
transforms and wavelet analysis routines subsequently 
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Fig. 2. — Zoom-in of a MBP group acting as the photospheric anchor point of a chromospheric spicule. Maximum longitudinal oscillatory 
power is found at the center of the MBP (green pixel; left panel), equating to power in excess of 7 X 10^ times the quiescent background. 
Phase shifts (in degrees), relative to this local maximum, are displayed using the red/blue color table, with maximum values in this example 
of +40 and —49 degrees. The right panel displays the phase information contoured on top of the simultaneous G-band intensity image, 
revealing how oscillatory power is closely confined within the MBP perimeter. 



applied on the resulting intensity fluctuations. We find 
significant oscillatory power (exceeding 500 times the 
quiescent background) in the averaged intensity of these 
G-band features, which typically occupy ~50 pixels, or 
130 000 km^. Furthermore, the most powerful intensity 
oscillations found in each photospheric MBP are at twice 
the period (130 - 440 s) of the corresponding chromo- 
spheric transverse oscillation, indicating the presence of 
frequency doubling. The dominant period remains con- 
stant over each MBP surface, indicating the presence of 
powerful, yet coherent drivers. We further utilize our 
high-resolution observations to deduce additional param- 
eters related to the longitudinal photospheric drivers, 
and find that the peak longitudinal oscillatory power is 
located at the center of each photospheric MBP, with a 
magnitude exceeding 7 x 10^ times the background av- 
erage (right panel of Fig. [2]). Power decreases radially 
away from the center of the photospheric MBP, eventu- 
ally dropping to zero at the edge of the structure. We 
compute the spatial distributions of oscillatory phase for 
the photospheric MBPs relative to this position of peak 
power. Figure [2] reveals how deviations in phase angle, 
from the positions of peak-power, are within the range 
of -|-40 to —49 degrees, providing absolute phase shifts of 
^90 degrees across each entire MBP structure. These 
phase shifts remain constant throughout the duration 
of periodic behaviour, with oscillatory power eventually 
ceasing as the MBP fragments. In addition, the detected 
photospheric MBP periodicities ( 130 - 440 s) overlap 
with the so lar p-mode spectrum ()Leighton et al.l [T962I : 
lUlrichl[r970( ). suggesting these global modes provide the 
underlying driving force required to generate spicule os- 
cillations. 

To investigate the temporal variations in wave be- 
haviour, we calculate the oscillatory power at each time 
step, and compare with simultaneous measurements of 
the same spicule observed in different wavelengths. Our 
analysis shows that the power of lower chromospheric 



(Ca II K) transverse oscillations peaks after the pho- 
tospheric (G-band) longitudinal periodicities, yet before 
the upper chromospheric (Ha) oscillations. Figure [3] re- 
veals how a discernible time lag exists between three dis- 
tinct layers of the lower solar atmosphere. This implies 
that the observed chromospheric spicule oscillations are 
a result of upwardly propagating magneto- acoustic wave 
modes, generated at the solar surface, which couple into 
transverse kink waves in a region encompassing the pho- 
tospheric/chromospheric boundary. An important trend 
to consider is how the power of detected oscillatory be- 
haviour varies as a function of atmospheric height. Power 
spectra for the G-band intensity, and Ca II K and Ha 
transverse oscillations, were normalised to their respec- 
tive quiescent backgrounds before comparison. Figure [3] 
shows an appreciable decrease in oscillatory power as a 
function of atmospheric height. This has important con- 
sequences, as it suggests the presence of physical damp- 
ing, and therefore conversion into heat, of the observed 
wave modes. Further examination of each MBP reveals 
multiple power peaks (also visible in Fig. [S]), suggest- 
ing the creation, and subsequent transverse behaviour of 
solar spicules, is quasi-periodic in nature. This is the 
first observational evidence of longitudinal-to-transverse 
mode coupling occurring in the Sun's atmosphere. 

To further investigate the propagation and cou- 
pling of these wave phenomena, we utilized advanced 
two-dimensional magneto-hydrodynamic (MHD) simu- 
lations. Thes e simulations wer e performed using the 
Lare2D code ()Arber et al.l |2001T ). where a spicule was 
modelled as a density enhancement, embedded in a 
straight, uniform magnetic field. The spicule density 
has a symmetric Epstein profile (Na kariakov fc RobertsI 
|1995[) in the transverse direction, with a peak density 100 
times the background value, and a diameter of approxi- 
mately 600 km (upper-left panel of Fig 14]). The magnetic 
field strength is chosen to give an external Alfven speed 
of 362 kms~^, and an Alfven speed of 36.2 kms~^ at the 
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Fig. 3. — Oscillatory power of an isolated MBP group in G-band 
(black), and associated Ca II K (blue), and Ha (red) spicule time 
series. The x-axis is time from the start of the observing sequence, 
while the y-axis is power averaged over the entire MI3P/spicule 
structure, and normalized to the background quiescent Sun. G- 
band power is for oscillations in intensity at a period of 258 s, while 
the Ca II K and Ha power is derived from the periodic transverse 
displacements of the spicule with a periodicity of 129 s. It is clear 
the G-band power leads the Ca II K and Ha time-series, indicat- 
ing the presence of an upward propagating magneto-acoustic wave, 
which couples into an transverse kink mode near the chromospheric 
boundary. 

centre of the spicule. A temperature profile is chosen to 
give a plasma (3 = 1, where eff icient mode coupling is 
expected to occur (IDe Moortel e t al. 2004). The simula- 
tion has a resolution of 600 x 2000 grid points, providing 
a numerical domain size of 4000 x 20000 km^. For the 
purposes of this investigation, only effects up to a few 
thousand km (i.e. a region incorporating the photosphere 
and chromosphere) will be considered. 

An artificial driver is applied at the lower boundary 
in the form of a longitudinal velocity component. This 
is consistent with the observational interpretation that 
the driver arises from a global p-mode oscillation. The 
coherent driver is applied across the diameter of the 
spicule, with a periodicity of 215 s, and a maximum am- 
plitude of 12.5 kms~^. However, to remain consistent 
with parameters derived from the observational time se- 
ries, one side of the simulated MBP has a phase dif- 
ference of 90 degrees with respect to the other (upper- 
left panel of Fig. 3]). As the driver is longitudinal and 
compressive, a 90 degree phase shift creates transverse 
(i.e. across the spicule) gradients in pressure which dis- 
place the spicule axis, generating the kink mode, in ad- 
dition to producing periodic compressions and rarefac- 
tions (sausage mode). These periodic pressure differ- 
ences, across the body of the spicule, also induce a fre- 
quency doubling of the coupled transverse wave, visible 
in the upper-right panel of Figure |4] by a change from 
a purely sinusoidal waveform. Both transverse kink and 
sausage modes are readily apparent in the lower panels of 
Figure m where time-distance cuts of observational, and 
simulated data, are displayed (see also Supplementary 
Movies 3 & 4). In our simulations, the magnitude of the 
excited kink mode is smaller than the amplitude of the 
input longitudinal driver. This is a consequence of the 
uniform magnetic field strength used. In the real solar at- 
mosphere, the magnetic field strength decreases rapidly 
with height, which is a likely contributor to the more 
prominent increase in kink amplitude seen in current 
ground- and space-based observations of chromospheric 





Fig. 4. — Upper left: The spicule density profile (solid line) nor- 
malised to the background value, while the dashed line shows a 
snapshot of the longitudinal velocity amplitude, 130 s into the sim- 
ulation, and scaled for clarity. Upper right: The transverse velocity 
signal, measured at the spicule axis, as a function of time for two 
heights. The dashed and solid lines represent the signal obtained at 
heights of 500 km and 2000 km above the photosphere, respectively. 
A periodic transverse velocity signal, which increases in amplitude 
away from the photosphere, reveals the progression of longitudinal- 
to-transverse mode coupling. Frequency doubling, represented by 
a change from a purely sinusoidal waveform, is observed in the sig- 
nal from a height of 2000 km. Middle: A time-distance cut of a 
chromospheric spicule observed in Ho. Bottom: A time-distance 
cut of the simulated spicule density, at a height of 2000 km above 
the photospheric boundary. There is a remarkable degree of sim- 
ilarity between the two lower panels, with both kink (transverse 
displacement of spicule axis) and sausage (periodic compression 
and rarefaction) modes visible. Supplementary Movies 3 & 4 dis- 
play time-lapse image sequences of the simulations presented here. 



spicules fe.g. iHe et al.l 120091 ). Nevertheless, the upper- 
right panel of Figure |4] indicates that the longitudinal- 
to-transverse mode-coupling mechanism present in our 
MHD simulations does promote a progressive increases 
in kink amplitude as a function of atmospheric height. 

By combining observed photospheric periodicities and 
phase distributions, in conjunction with realistic MBP 
structuring, our simulations reveal efficient, and ubiq- 
uitous, mode conversion of longitudinal oscillations into 
their transverse kink counterparts, at twice the driven 
frequency. The simulations display wave properties con- 
sistent with our observations, and demonstrate how 
transverse spicule oscillations, with significant ampli- 
tudes, readily exist in the solar atmosphere, and there- 
fore have important consequences for energy transporta- 
tion into the solar corona. We estimate the energy 
flux, E, of the observe d chromospheric waves using 
()De Pontieu et al.|[2007b[) . 



'El — pv VA 



(1) 



where p is the mass density of the flux tube, v is 
the observed velocity amplitude and va is the corre- 
sponding Alfven speed, defined as va = B/^fJi^p, with 
Pq the magnetic permeability. For a mass density of 
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p PS 1.3 X 10^* kgm"'^, and a magnetic field strength 
B Rs 12 G, derived from a br ight network chromospheric 
model ([Vernazza et al.iri981f) . coupled with an observed 
velocity amplitude of f « 15 kms^^, the energy flux 
in the chromosphere is i? w 3 x 10^ Wm^^. It is es- 
timated that appr oximately 2.2% of the solar su rface is 
covered by MBPs (|Sanchez Almeida et ani201Clf ). and if 
each MBP is linked to a single corresponding chromo- 
spheric Type I spicule, it equates to an average global 
energy of 660 Wm~^. Current work suggests waves 
with an energy flux «100 Wm~^ are sufficiently vig- 
orous to heat the localized corona and launch the sola r 
wind when their energy is thermalized (|Jess eTaDHooi). 
Therefore, a transmission coefficient of «15% through 
the thin transition region would provide sufficient energy 
to heat the entire corona. Regions on the solar surface 
containing highly magnetic structures, such as sunspots, 
pores and large MBP groups, should possess even higher 
mass densities and magnetic field strengths, in addition 
to a greater numbers of spicules. In this regime, the en- 
ergy fiux available to heat the corona will be significantly 
higher than the minimum value required to sustain local- 
ized heating. 



4. CONCLUDING REMARKS 

We have utilized images of high spatial and tempo- 
ral resolution, obtained with the Rapid Oscillations in 
the Solar Atmosphere (ROSA) instrument at the Dunn 
Solar Telescope, to reveal how the generation of trans- 
verse kink oscillations in Type I spicules is a direct re- 
sult of mode conversion in the lower solar atmosphere. 
Through comparison of our observations with advanced 
two-dimensional magneto-hydrodynamic (MHD) simula- 



tions, we show how longitudinal pressure modes in pho- 
tospheric magnetic bright points (MBPs), with period- 
icities in the range 130 - 440 s, funnel upwards through 
the Sun's atmosphere, before converting into kink modes 
at twice the initial frequency, often with amplitudes ex- 
ceeding 400 km. We conclude that the mode conver- 
sion and period modification is a direct consequence of 
the 90 degree phase shift encompassing opposite sides of 
the photospheric driver. This is the first observational 
evidence of this mechanism occurring in the solar at- 
mosphere. Energy flux estimates for these oscillations 
indicate that the waves are sufficiently energetic to ac- 
celerate the solar wind and heat the quiet corona to its 
multi-million degree temperatures. 

The naming of transverse oscillations observed in so- 
lar structures remains highly controversial, with def- 
initions revolving around "Alfven", "Alfvenic", and 
"magneto-sonic kink" terminology. While our observa- 
tions demonstrate signatures consistent with previous 
studies on Alfvenic waves (e.g. iDe Pontieu etal]|2007bl: 
iMcIntosh et al.|[2bllD , we have deliberately chosen to de- 
scribe the observed periodic motions simply as transverse 
kink waves. This is currently the most unopposed de- 
scription of such wave phenomena, and avoids the poten- 
tial pitfalls of what is a rapidly developing area within 
solar physics. 
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